EVALUATION OF DEDEKIND SUMS, EISENSTEIN COCYCLES, 
AND SPECIAL VALUES OF L-FUNCTIONS 



PAUL E. GUNNELLS AND ROBERT SCZECH 

Abstract. We define certain higher-dimensional Dedekind sums that generalize 
the classical Dedekind-Rademacher sums, and show how to compute them effec- 
tively using a generalization of the continued-fraction algorithm. 

We present two applications. First, we show how to express special values of 
partial zeta functions associated to totally real number fields in terms of these 
sums via the Eisenstein cocycle introduced by the second author. Hence we obtain 
a polynomial-time algorithm for computing these special values. Second, we show 
how to use our techniques to compute certain special values of the Witten zeta- 
function, and compute some explicit examples. 



1. Introduction 

1.1. Let a he a square matrix with integral columns cij G (j = 1, . . . , n), and let 
L C Z" be a lattice of rank > 1. Let v G M", and let e G with Cj > 1. Associated 
to the data [L, a, e, v) is the Dedekind sum 

(1) ^ = ^(L,a,e,t;) := V e((x,t;))- r^^, r^- 

Here (x, y) := ^ Xii/i is the usual scalar product on M", e(t) is the character exp(27rit), 
and the prime next to the summation means to omit terms for which the denominator 
vanishes. This series converges absolutely if all Cj > 1, but only conditionally if cj = 1 
for some j. In the latter case, we define the value of S by the Q-limit 

(2) a{x) ^■.= \im(^ ^ a{x] 

\Qix)\<t 

where Q is any finite product of real-valued linear forms on that do not vanish 
on Q" \ {0}. As we explain in § p^.6| , this limit depends on Q in a rather simple way. 
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-2TTi) 



Nevertheless, to keep notation to a minimum, we assume for now that the series (|T]) 
converges absolutely. 

1.2. The arithmetic nature of the values of S is well known. Up to a power of 
27rz, they are always rational numbers if v e Q". However, the explicit calculation 
of these values is not easy, especially if | deto"! is large. In this paper, we exhibit 
a polynomial-time algorithm for calculating S efficiently. Before stating our main 
result, we briefly review a few special cases. 

In the case n = 1, a = 1, L = Z, we have 

(3) 5(L,a,e,t;) = -^^^,(z;), 

e! 

where ^eiv) is the periodic function (with period lattice Z) that coincides with the 
classical Bernoulli polynomial Be{v) on the interval < w < 1. 
More generally, if L = Z", then we show in Proposition that 

(4) S{L,a,e,v) = K ^ ^^^iui) ■ ■ ■ ^eA^n), 

reL/aL 

where 

n 

(5) u = a^^{r + v) and k = sgn(det a) JJ^ 

The flnite sum on the right of (^) is a classical Dedekind sum. Although theoretically 
satisfying, from a computational point of view this sum is of interest only if | deto"! 
(the number of summands) is relatively small. 

1.3. In general, if rankL < n, no simple analogue of (^ seems to exist except for 
the special case where rankL equals the number of linear forms {x,crj) in (|I]) that 
are not proportional to each other when restricted to the subspace L M. Dedekind 
sums of this special type are called diagonal. It is clear that the flnite sum formula 
@) remains valid (modulo obvious modiflcations) for all diagonal Dedekind sums. 

Let us further say that a diagonal Dedekind sum is unimodular if the corresponding 
flnite sum has exactly one term. As a measure of the deviation from unimodularity, we 
introduce an integer- valued function || ||, called the index, on the set of all Dedekind 
sums (Deflnition p.6|) . For example, if L = Z", then \\S\\ = \ deta\. In particular, S 
is said to be unimodular if and only if 1 15*11 = 1. 

1.4. We are now ready to describe our main result. 

1.5. Theorem. Every Dedekind sum S{L, a, e, v) can be expressed as a finite rational 
linear combination of unimodular diagonal sums. If n, rankL, and e are fixed, then 
this expression can be computed in time polynomial in log US'!! . Moreover, the number 
of terms in this expression is bounded by a polynomial in log \\S\\. 
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The special case n = 2 and e = (1, 1) has been known for a very long time. This 
is the case of the classical Dedekind-Rademacher sums, where the simple form of 
their reciprocity law combined with the Euclidean algorithm immediately yields a 
polynomial-time algorithm for computing their values (cf. Example [7.5| ). Our proof 
of Theorem p. .51 is similar. First we prove a general reciprocity law for the sums 
S. Applying and specializing this law repeatedly yields a representation of 5* as a 
linear combination of diagonal sums. Combining it further with the algorithm of 
Ash-Rudolph [0 finally yields a representation by unimodular diagonal sums. 

1.6. Among the many applications, we review in this paper the problem of calcu- 
lating the Eisenstein cocycle [0] on GL„(Q) that, when combined with Theorem |L5 



yields a polynomial-time algorithm for calculating special values of partial zeta func- 
tions of totally real number fields (§§§]-0). These special values are of great interest in 
view of the many conjectures they are subject to (Leopoldt conjecture, Brumer-Stark 
conjecture, etc.). We also review the connection between Witten's zeta function and 
Dedekind sums, and use our techniques to give some explicit formulas for special 
values of these functions (§|^). 

1.7. Acknowledgements. We thank Gautam Chinta for carefully reading prelimi- 
nary versions of this article, and for help with the examples in §^ 

2. The modular symbol algorithm and Dedekind reciprocity 

2.1. Let (Ti, . . . , (T„ G Z" be nonzero primitive points, and \et D = \ det((Ti, . . . , ct„)|. 
Let tf G Z" be another nonzero primitive point, and define 

Di{w) := I det(ai, . . . ,ai, . . . ,ct„,w)|, i = l,...,n. 

The following basic result will play a key role: 

2.2. Proposition. [|1], If D > 1, then there exists w E ^ {0} such that 
(6) < D,{w) < i = l,...,n, 

and at least one Di{w) ^ 0. Moreover, for fixed n the point w can be constructed in 
polynomial time. 

Proof. (Sketch) Here we only show that w exists. Let P be the open parallelotope 

Then P is an n-dimensional centrally symmetric convex body with volume 2". By 
Minkowski's theorem (cf. [0, IV.2.6]), P fl Z" contains a nonzero point. This is the 
desired point w. 

□ 
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2.3. Remark. Ash-Rudolph |1|] show that w satisfies < Di{w) < D. They also show 
how to construct w using the Euclidean algorithm. The stronger estimate (^) and the 
statement about polynomial time are due to Barvinok |Q, Lemma 5.2]. In practice, 
to efficiently construct w such that Di{w) is small, one may use LLL-reduction |]^, 
§3]. 

2.4. We now state the Dedekind reciprocity law. For any nonzero point v G M", let 
be the hyperplane {x \ {v,x) = 0}. Let Q be a finite product of real-valued linear 

forms on that do not vanish on Q" \ {0} and recall that the notation 



a[x 



means the sum is to be evaluated using the Q-limit (§|1.1|, §pT^). 

2.5. Proposition. Let ctq, . . . , o"„ G be nonzero. For j = 0, . . . , n, let be the 
matrix with columns ctq, . . . , a^, . . . , cr„. Fix L C Z", and assume e = 1 := (1, . . . , 1). 
Then for any v G M", we have the following identity among Dedekind sums: 



j=0 



(7) 

Proof. Let Dj = det . We have an identity for rational functions of x 



j=0 



n 



7=0 nfc^j' 



0, 



valid for any x G satisfying {x, (Tj) 7^ for j = 0, . 
(n + 1) X (n + 1) matrix 



, n. To see this, consider the 



(9) 



(x, cTo) ... (x, a„) 

CTq ... (T„ 



This matrix is singular, since the first row is a linear combination of the others. 
Expanding by minors along the top row, and dividing by Yli^y ^k) yields (|^). 

To pass from (|) to (^, we need to incorporate the exponential character and sum 
over L using Q. There is no obstruction to doing this, although we must omit terms 
where any linear form vanishes. We obtain the expression 



(10) 



5:(-lK^e({x,.))^ 



0, 



Q 



where the inner sum is taken over all x G L with (x, ak) 7^ for A; = 0, . . . , n. 

In (|10|), the jth inner sum corresponds with the Dedekind sum S{L,a^, 1) except 
for the terms with (x, aj) = and (x, cr^) 7^ for k 7^ j. In other words, to make the 
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jth sum into a Dedekind sum, we must add 

Simultaneously adding and subtracting (pA]) to (||) yields □ 
2.6. We now recall the Q-limit formula from [0, Theorem 2]. Let 

m 
i=l 

be a product of m > 1 linear forms 

n 

Qiiy) = ^Qijyj^ 

with rationally independent real coefficients Qij. We think of Q as an m x n matrix 
with rows Qi. 

Given a vector e = (ei, . . . , e„) of positive integers and a vector v G M", let 

J = {j I = 1 and Vj G Z}. 

If # J = mod 2, define 



(12) B.(.,Q)=i-f:(n^)n^..(..) 

otherwise let Be(f , Q) = 0. In particular, if J = 0, then 



i=i 

Now the Q-limit formula can be stated as follows. Let Id„ be the n x n identity 
matrix. Then 

(13) 5(Z^Id„,e,^;) =kM,{v,Q), 

Q 

where 



(14) K = n 



■ 1 
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3. DiAGONALITY AND UNIMODULARITY 

3.1. We begin with some simphfying assumptions to ease the exposition. We define 
the rank of S = S{L, a, e, v) to be the rank of the lattice L. 

Suppose that the rank of S is i, and for any k let Z'' C Z" be the sublattice 
spanned by the first k standard basis vectors. We claim that S can be computed 
using Z^ instead of L. Indeed, writing L = gZ^ with a matrix g G GL„(Q), and 
letting cr' = g^a, v' = g^v, and Q' = Qg, we have 

S{L, 0", e, v) = (det g)^^S{Z^, a', e, v') 
Q 

The entries of a' need not be integral, but after multiplying by an appropriate rational 
factor we can assume this is true. In fact, by further multiplication of a' by rational 
numbers and permuting columns, we can write 



S{L, a, e, v) = qS{Z^, a', e, v') 



q e 



where the pair {Z^, a') satisfies the following: 

(i) For each column cr^-, the vector a'^ fl Z^ is primitive and integral. 

(ii) If two columns of a' induce proportional linear forms on Z^, then these two 
linear forms coincide on Z^, and are adjacent columns of a'. 



3.2. Definition. We say that a rank i Dedekind sum is normalized if the conditions 
above are met. 

3.3. Let S{Z^, a, e, v) be a normalized Dedekind sum, and let = ^ ej. We claim 
that without loss of generality, we need only consider sums for which e = 1. Indeed, 
let i?" M.^ be the spans of the first £ (respectively n) basis vectors in M^, 
and identify M" with R^. Let Q' be a product of linear forms on such that Q' 
restricted to R^ equals Q restricted to M". We claim that we can construct an x 
matrix a' and v' G such that 



(15) 



S{Z\a,e,v) = S{Z^,a',l,v') 



To see this, let vr: — > be the projection on the first i components. Given 
V G -R", call any v G such that n^v) = 7r(f ) a lift of v. 

Now to construct S{Z^,a', l,v'), let v' be any lift of v. If the restriction of the 



linear form aj to Z appears on the left of (|T^) with multiplicity ej, set j columns of 
a' to be j different lifts of aj. If we choose these lifts so that det a' = det cr, then we 
obtain (1151). 



3.4. Definition. We say that a Dedekind sum with e = 1 is properly embedded. 
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3.5. Let S{Z^, a, l,v) be a properly embedded, normalized Dedekind sum. Let |n] 
be the set {1, . . . , n}. 

3.6. Definition. The index of S, denoted H^H, is defined to be 

Max|det(7r((7iJ,...,7r((7iJ)|, 

/cH 

where the maximum is taken over all subsets / = {ii, . . . ,ii} of cardinality i. A 
Dedekind sum is called unimodular if ||5'|| = 1. 

3.7. Now define a partition 

s 

(16) M = U^fc' ^<^<^ 

fe=i 

as follows. Put 

i, J G Ik if and only if 7r((Tj) = 7r((Tj). 

In other words, two elements of \n\ are in the same set of the partition if the corre- 
sponding columns of a induce the same linear form on Z^. 

Let Pk — H^Ik- The vector p{S) = {pi, ■ ■ ■ ,Ps) is called the type of S. To emphasize 
the type, we relabel the columns of a as 

(17) ial,...,al\al,...,a',\...,al,...,a^/). 

For any k — 1, . . . ,s and any i — 1, . . . ,pk, we denote the point 7r((j^) e by o"/^. 

3.8. Definition. A Dedekind sum is called diagonal if p{S) has length i. 

We omit the proof of the following simple lemma: 

3.9. Lemma. Any normalized, properly embedded rank 1 Dedekind sum is both di- 
agonal and unimodular. 

3.10. The virtue of diagonality and unimodularity is the following: 

3.11. Proposition. Keep the notation of the preceding section. Let S{Z^,a,l,v) be 
properly embedded, normalized, and diagonal. Let p be the £ x i matrix (u/j, . . . , cr/J, 
and let (pi, . . . ,pi) be the type of S. Then S\q is well-defined. Moreover, 

(18) SiZ',a,l,v) 

Q detp 

where 

(19) M = p"^(r + 7r(t;)), = Q" o p-* o tt, and 

-27ri)P^ 



(20) 
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Here Q" is the restriction of Q to into . If S is unimodular, the sum (|T8|) has 
only one term. 



Proof. The second statement follows easily from the first, so we focus on the first. 
By definition, we have 

*3 ^^ryp 111 



(21) 



ll<k<i 



{x,al) 



Pk 



Letting Q' = Q" o p * o vr and v' = p ^7r(f), the right of (l2l|) becomes 
(22) detp 



det(a) $^'e((y,t;'))^ 



Pk 



Inserting the character relations 



(23) 

we obtain 
(24) 

(25) 



^ e((?/,p V)) 



det 



a 



I det p\ 

det cr 
I det p\ 



E E' 



e((?/,p ^(r + t;))) 



5^ s(z^Id.,p, 



where u = p ^(r + vr(t')). 

The proposition follows now from the Q-limit formula (^ 

4. Algorithms 



□ 



4.1. In this section we prove that any Dedekind sum is a Q-linear combination of 
diagonal, unimodular sums. We begin with a lemma. For simplicity, we abbreviate 
the Dedekind sum to S{Z^,a). 

4.2. Lemma. Let S{Z^, a) be a normalized, properly embedded Dedekind sum, and 
let cTj be a column of a. Then 

\\S{Z'na^,a)\\<\\S{Z\a)\\. 
Proof. Without loss of generality, assume that (Xj = (Xi. We can represent a as 



(26) 



* 



0-/2 



where there are pk columns of the form {(Ti^, *)*. The stars in the last n — i rows 
represent numbers that are irrelevant, since they don't affect the value of the sum. 
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Let 7 G GL„(Q) be a matrix that carries fl o\ onto ^ . For = 1, 



let 



7cr^, where a\ is any hft of Oi^,. Then 70" has the form 



(27) 



70" 








■■■ 




■■■ a\ 







... 




... 


* 


* 


* 




* 




1, 1, 


and n — £ 


and 


e = ±1. 



Here the row blocks have sizes 

Now let d be the determinant of any (i — l) x (i — l) minor from the top i—1 rows 
of p7|). This determinant will be the same up to sign as the determinant of an i x i 



minor containing (Ji from the top i rows of (pGl) . Hence \d\ < \\S\\, and the proof is 
complete. 



□ 



Now we come to our first main result. 



4.3. Theorem. Let S = S{Z^, a) be a normalized, properly embedded Dedekind sum. 
Then we may write 

(28) S = J2 Q,S{Z', g) + J2 qrS{Z'-\ r), 

Q&R reT 

where 

1. the sets R and T are finite, 

2. q,, qr e Q, 

3. each S{Z^, g) is diagonal, and 

4. each of the Dedekind sums on the right of (f 



Proof. Write a = {al, . . . ,a{^,a. 
columns if necessary, we may assume that u/^ , . . . , aj^ are linearly independent in 
Z^. We will show that we can write an expression like ( ^8]) so that the rank i sums 
on the right have type 

{pi,...,Pi-l,...,pe,...,Ps + l), 

for some 1 < i < i. Iterating this construction proves that we can have the rank i 
sums on the right of (^8]) diagonal. 

We proceed as follows. Since ct/^ , . . . , o"/^ are linearly independent, we can find 
unique rational numbers aj such that 



^P2 



has index < \\S\\. 
■ jCTg") as in (|l^). Permuting 



(29) 



Now we want to apply the relation in Proposition |2.5| . For each o"/^. , we choose pj 



rational lifts aj, 



not necessarily equal to the columns of a, and use them to 



form a matrix a. Clearly S{Z ,a) = S{Z ,a). 
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Now define 

Clearly 7r(w) = cr/^. Write al {w) for the matrix made from a by replacing the column 
al with w. Using the columns of a and w m. we find 

(30) S{Z', a) = S{Z' n w^, a) + J2(^l)S{Z', aHw)) 

where el G {±1} is determined by the reciprocity law, and the sums are over pairs 
(i, j) satisfying 1 <i < s and 1 < j < Pi- 

We claim that the sums in (^) actually have only ^ terms. This follows since the 
points 

are dependent. Hence any sum such that al (w) contains these points vanishes. More- 
over, the sum S{Z^ fl vj-^, a) is zero, since by construction a column of a induces a 
linear form vanishing on Z^ (Iw^. Hence ( ^O]) becomes 

e e 

(31) S{z',a) = Y,s{z',al{w)) -Y,s{z'n{al)\al{w)). 

i=l 1=1 

Now consider the types of the rank i Dedekind sums on the right of (^T]). If the 
type of 5* was 

then the type of S{Z^,a}{w)) is 

{pi,...,Pi-l,...,Pe,...,Ps + l). 

Hence by induction we can write S" as a finite Q-linear combination of diagonal rank 
i Dedekind sums plus sums of lower rank, which proves 1-3 of the statement. 

To complete the proof, we must show that the indices of the Dedekind sums on 
the right of (|28| ) are no larger than \\S\\. Indeed, Lemma ^]2| implies that the indices 
on the right of (|30|) are no larger than \\S\\, and so the claim follows. □ 

4.4. Remark. There are some simplifications in ( pT]) that are worth mentioning if one 
wishes to implement the diagonalization algorithm. First, if we construct a so that 
det 5" = 1, then deta^?^ = aj. This means that these determinants can be computed 
when one computes aj^. Second, not all the terms on the right of (|3lD necessarily 
appear. In particular, the rank i — 1 sum S(^Z^ fl {crl)^, al{w)) appears on the right 
of (0) only iipi = l. 
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4.5. Theorem. With the notation as in Theorem we can write 
(32) S = J2 %S{Z', q) + Y, qrS{Z'-\ r), 

g£R reT 

where the sums of rank i on the right are diagonal and unimodular, and the rank i — 1 
sums have index < \\S\\. 



Proof. By Theorem |4.3|, we may take S to be diagonal. Suppose that 

D = I det((j/i, . ..,aij\ > 1. 



Then by Proposition p.2| , there exists w & such that 

Dj = I det(a/i , . . . ,aj^, . . . ,ai^,w) 



satisfies < Dj < D^^ As in the proof of Theorem write w = Ylj=i ^j^iji 
and set w = Yl^j=i ^j'^j- 



Now apply Proposition |2.5| , using w and the columns of a, to write S* as a finite 
Q-linear combination of new Dedekind sums. These sums won't be diagonal, but we 
can apply the proof of Theorem [4.3| with w playing the role of aj^ . The resulting sums 
will include lifts of w in their columns and will be diagonal. Thus the resulting rank 
i sums will have index < \\S\\, and the sums of lower rank will satisfy the conditions 
in the statement of Theorem [4.3| . By induction on the index, this completes the 
proof. □ 

4.6. Corollary. Any Dedekind sum can be written as a finite <Q-linear combination 
of diagonal, unimodular sums. 



Proof. First normalize and embed properly. By Lemma |3]^, any rank 1 Dedekind sum 
is automatically unimodular and diagonal. The result follows by applying Theorems 
[4.3| and [4.5| and descending induction on the rank. □ 

5. Complexity 

5.1. In this section we discuss the computational complexity of Corollary [4.6| . In 
particular, we show that if n and £ are fixed, and 5* = S{Z^, a) is normalized and 
properly embedded, then we can form a finite Q-linear combination of diagonal, 
unimodular Dedekind sums 



S = J2<l,S{Z\g), 



where is bounded by a polynomial in log \\S\\. As a corollary we obtain that this 
expression can be computed in polynomial time. 

To do this, we must make a more detailed analysis of proofs in We begin by 
analyzing diagonality. 
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5.2. Lemma. Let S = S{Z^, a) be a normalized, properly embedded Dedekind sum, 
and write 

(33) S = J2 %S{Z', q) + Y. qrS{Z'-\ r) 

as in Theorem so that the rank I sums in (|33| ) are diagonal. If i > 1, there exist 
constants Mn/ and Nn/ such that 

i^R<Mn,e and #r < iV„,,. 

Proof. Write p{S) = {pi, . . . ,Ps), where i < s < n. By the proof of Theorem O, we 
know how to pass from a sum of type 

(34) {pi,...,pi,...,pe,...,ps) 
to a hnear combination of sums of types 

(35) {pi,...,pi-l,...,pi,...,ps + l), i = l,...,i. 

By iterating this, we pass from S* to a hnear combination of sums with types 

(36) {p[,...,p[_^,0,p[_^^,...,p[,...,p'^), i = l,...J. 

We win bound the number of rank £ (respectively rank i — 1) sums produced in 
passing from (|3^) to (|36D by a constant M^j (resp. N^j). We can then take 



(s) 

i+1 



and similarly for Nn,£- 

To describe what happens in going from (^) to (^), we use a geometric construc- 
tion. Let B = B{pi, . . . ,pi) he the set 

B= {{xu...,xe) el} I Q<Xi<pui = l,...,l]. 

The points in B correspond to types of intermediate sums in the passage from (B^) 
to (|36D . In particular, passing from (|3^) to (^) can be encoded by moving from 
(xi, . . . , x^) to (xi, . . . , Xj — 1, . . . , x^) in B. Moreover, the sums of the form (|36D 
correspond to the subset of points Bq C B with exactly one coordinate 0. (See 
Figure |l].) 

Now the constant M^^j will be given by Max#i?o, as i?(pi, . . . ,pe) ranges over all 
possibilities for fixed n, i, and s. If we allow the pi to become continuous parameters, 
then a simple computation shows that the maximum occurs when pi = ■ ■ ■ = pe = 
(n — s + With these conditions we have 
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Figure 1. The set 5(2,3). The arrows represent apphcations of 
Proposition ^ The black dots represent points in Bq. 

The constant A^^*] can be computed similarly. There are (-^ + 1) sums of rank 
produced for each point in 

B+ := {{xi,...,xe) eB\xi^O}. 

One finds again that the maximum occurs when pi = ■ ■ ■ = pi = {n — s + 



n — s + 



and is 



□ 



5.3. Proposition. // n and i are fixed, then ( |33D can be constructed in constant 
time, independent of \\S\\. 

Proof. This follows easily from the proof of Lemma |5.2| . Forming the expression (|33[) 
is purely combinatorial, and makes no reference to \\S\\. In particular, the number of 
steps needed can be bounded for fixed n and i. □ 



5.4. Now we investigate the size of the output in Theorem [4.5| . The main step of 
the proof of Theorem shows how given S, one may write 

(38) S = Y, <1,S{Z', 0) + Y. <lrS{Z'-\ r), 

where the sums of rank i on the right satisfy \S[Z\ q)\ < Denote by d 

(respectively Ce~i) the number of rank i (resp. rank i — sums on the right of (p8|). 



5.5. Lemma. 



We have Cf, < M^^+^^J and < A^^^+^^ 



n+l,. 



Proof. The proof is very similar to the that of Lemma |5.2| , with the following twist: 
we begin with a diagonal sum, increase the number of distinct linear forms by one, 
and then make the sums diagonal again. We can keep track of the number of sums 
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produced using the set B as in the preceding proof, although we must replace n with 
77, + 1 to accommodate the extra initial step. We leave the details to the reader. □ 



5.6. Now consider the expression (|38D. To complete the proof of Theorem [4.5| , 
we repeat the process that produced ( |55D until all rank i sums are diagonal and 
unimodular, and we obtain the expression (^). Using Lemma ^ and the estimate 
(|^), we can bound the number of rank i sums produced. 

5.7. Proposition, (cf. 0, Theorem 5A]) Write S as a sum of diagonal, unimodular 
rank i sums and lower rank sums as in (|32D . Then the number of rank i sums on the 
right of (^) is bounded by 

(39) C'(log||S||)^l°^^^ 

where A = \og{n/{n — 1))~^, and C is a constant independent of S. 

Proof. In (^SD, we have 

Q<\\S{Z\q)\\<\\S\\^--'^I\ 
Thus after t iterations we'll have 

< \\S{Z\q)\\ < 

Since the index of a Dedekind sum is always an integer, the condition for termination 
is that for some £ > 0, we have 

loglogll^ll -loglog(2-£) 



(40) II 5|| (("-!)/«) < 2-e, or t> 



log n — \og{n — 1) 



On the other hand, by Lemma |5.5| , we know that t iterations will produce no more 
than C\ sums of rank £. So fix e > 0, set A = logn — log(n — 1), and let 

^, //-loglog(2-e) \, ^ 

c; = exp M — — + ij log a 

Then 

,A\ogCi_ 



C\ < C',{\og\\S\ 

Now if we define C = Max£<„C^, we obtain (|39|). □ 

5.8. We are now ready to discuss the complexity of our algorithms. 

5.9. Theorem. Let S = S{Z^, a) be a normalized, properly embedded Dedekind sum. 
Using Corollary \^.(\ , write S as a Q-linear combination of diagonal unimodular sums. 
Then there exists a polynomial Pn/ such that the number of terms in the expression 
is bounded by Pn/{\og\\S\\) . Moreover, we have 



deg Pn/ < A log 
where A = log{n/[n — l))^"*^. 



2132... £(^-1) 



EVALUATION OF DEDEKIND SUMS 



15 



Proof. Fix n. We proceed by induction on i. 

First, if £ = 1, then by Lemma |3.9| the sum S{Z^,a) is aheady diagonal and 
unimodular. Hence we may take P„ i = 1. 

Next, assume that the statement is true for sums of rank i — 1, and let Pn,e-i be 
the corresponding polynomial. First we claim that without loss of generality, we need 
only consider the case that S is diagonal. Indeed, apply Theorem and write 



(41) 



£-1 



geR 



reT 



where the rank i Dedekind sums are diagonal, and all the Dedekind sums have 
index < \\S\\. By Lemma f).2[ the sets R and T have a bounded number of elements 
independent of 5*. Hence we may bound the output for diagonal S, and then multiply 
this by a constant to obtain our final answer. 

Now we apply Theorem |4.5| , and we must count the number of rank Dedekind sums 
produced. By Proposition |5.7| , we know that the total number of rank i sums will be 
bounded by 



(42) 



C'(log 



AlogCe 



Furthermore, each sum of lower rank produced in the proof of Theorem can be 
written as a sum of < P„^^_i( 115*11) Dedekind sums by induction. To find the total 
output, we must count these lower rank sums. 

We can do this as follows. Let Q = Qn,e '■= C'^-i-Pn.^-i, where Ce-i is the constant 
in Lemma ^.5| . Represent the process of reducing the diagonal sum S to unimodularity 
by the following diagram: 



(43) 



rank 



rank i — 1: 




The top row represents the bound on the number of rank £ sums at each step of the 
algorithm, and the bottom row is the number of rank i — 1 sums. 
According to (H3), we have produced 



AlogCi 



a-1 



sums of ranks < i. Adding this estimate to (^), we find that the total number of 
Dedekind sums produced will be a polynomial of degree 



degQ + A log Ci = deg P„,^_i + A log 
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To complete the proof, we compute the degree of Pn/ by induction. Indeed, using 
the estimate 

and that degP„^i = 0, an easy computation shows that 

degP„,, < A log l^-^^—^ j , 
as required. □ 
5.10. Example. Here is a table of the bound of deg Pn/ for small values of n and £. 





2 


3 


4 


5 


6 


7 8 


2 


1 












3 


2 


5 










4 


4 


10 


15 








5 


7 


16 


25 


33 






6 


9 


23 


37 


50 


60 




7 


12 


30 


50 


69 


86 


99 


8 


15 


38 


64 


90 


114 


135 150 



5.11. Corollary. Keeping the same notation as in Theorem \5. % for fixed n and i 
we may express S as a finite Q-linear combination of diagonal unimodular sums in 
time polynomial in log \\S\\. 



Proof. First, the vector w constructed in Proposition p.2| can be found in polynomial 
time in the size of the coefficients of S. In fact, investigation of |Q, Lemma 5.2] 
shows that the rational numbers aj in (|29| ) in the proof of Theorem [4.3| can also be 



constructed in polynomial time. This implies that w and the can be found in time 
polynomial in log || 5*11 . The proof of Theorem then shows that the final expression 



can be computed in polynomial time. □ 



6. The Eisenstein cocycle 

6.1. In this section, we briefiy review the construction of the Eisenstein cocycle 
introduced in |^. In particular, we show that this is a finite object that can be 
calculated effectively using Corollary [4.6| . Roughly speaking, the Eisenstein cocycle 
represents a generalization of the classical Bernoulli polynomial within the arithmetic 
of the unimodular group V = G'L„(Z). 
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6.2. Let = {Ai, . . . , An) be an n-tuple of matrices Ai G GLn{^). For an n-tuple 
d = {di, . . . , dn) of integers 1 < (ij < n, let £^ (d) C M" be the subspacc generated by 
all columns A^j such that j < di. (Here Aij denotes the jth column of Ai.) Writing 
^(d)-*- for the orthogonal complement oi ^{d) in R", we let 

n 

(44) X{d) = ^{d)^ \ U a^, where ai = Ad,- 

i=l 

The n-tuple determines then the stratification 

(45) M"x{0}= □X(c^), 

indexed by the finite set 

D = D{^) = {d I X{d) ^ 0}. 
Associated to this decomposition is the rational function ip{s^) on M" \ {0} defined 

by 

^ det(cTi, . . . , a„) 

6.3. More generally, if P{Xi, . . . , X„) is any homogeneous polynomial, we form the 
differential operator P{—dxi, ■ ■ ■ , —dx^) in the partial derivatives dx^ :— d/dxi, and 
define 

^{S^){P, X) = P{-dx„ -dxMs^){x). 

The last expression can be written more explicitly as 

n 

(46) ^{^){P, x) = det(a) Pr{^) H T^^aIT^' 

where r runs over all decompositions of deg(P) = ri + • • • + r„ into nonnegative parts 
Tj > 0, and Pri^) is the homogeneous polynomial in the (7^ defined by the expansion 

p(x.*)=^p»n5^. 

r j=l 

In the excluded = 0, it is convenient to set iIj{s2/){P, 0) = 0. 

The definition of the Eisenstcin cocycle \E' is now easy to state: 

(47) <i/{^){P,Q,v) := (27ri)-"-^"s-P ^ e{{x,v))iP{^){P,x) 

The series on right converges provided all components Ai of ^ are in GL„(Q). 
However, since the convergence is only conditional, we are forced to introduce the 
additional parameter Q specifying the limiting process. 
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6.4. Let M be the set of all complex-valued functions f{P,Q,v) with P, Q, v as 
above {v G M"). Then M is a left F-module under the action 

Af{P, Q, v) = det{A)f{A'P, A-^Q, A-\), AeT, 

where the implied F-action on homogeneous polynomials is given by {AP){X) = 
P{XA). With respect to this action, the map \E': F" — > M has the property 

(48) ^(A^/) = Am{jz/), A e F, ^ e F'^, 

n 

(49) *(^0' • • • , • • • , A„) = 0, A, G F. 

In other words, \E' is a homogeneous cocycle on F. It is known that \& represents a 
nontrivial cohomology class in if"'~^(F; M) |]^, Theorem 4]. 

Combining (0)-(47), we see that is a finite linear combination of Dedekind 
sums, 

vI/(^)(P,Q,t;) = (27rO-"-'^<=^'^5^5^P,(a)5(L,a,e,' 

d&D r 

Here a is the matrix with columns Aj^- for i = 1, . . . , n, L is the lattice ^(rf)-*" fl 7/' 
and = 1 + Tj. The case P = 1 is of special interest: 

vl/(^)(l, g, t;) = {2m)-- ^(^' ^' 1' ^ 

del) 

This case yields the classical Dedekind-Rademacher sums if n = 2, and, more impor- 
tantly, it corresponds to special values of partial zeta functions at s = 0. 



7. Values of partial zeta functions 

7.1. Let P be a totally real number field of degree n over Q, and let f, b be two 
relatively prime ideals in the ring of integers ffp. The partial zeta function to the 
ray class b mod f is defined by 

C(b,f,s):= ^(^)"'' 3^(^)>1> 

o=b mod f 

where a runs over all all integral ideals m. ffp such that the fractional ideal ab~^ 
is a principal ideal generated by a totally positive number in the coset 1 + fb~^. 
According to Klingen-Siegel, the special values C,{b, f, 1 — s), where s = 1, 2, 3, ... , 
are well-defined rational numbers. In this section, we give a formula for calculating 
these numbers in terms of the Eisenstein cocycle \E'. 
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7.2. The formula depends on the choice of a Z-basis W for the fractional ideal 
fb"^ = ^Ziy,-, together with the dual basis W* determined by TT{W*Wj) = 6ij. 
Here we identify a G F with the row vector (a^^\ . . . , a^'^^) G M", where the a^-'^ 
are the n different embeddings of a into the field of real numbers. Then W can be 
identified with a matrix in GL„(R) whose jth row is the basis vector Wj. Let 

P{X) = N{b)l[J2x,wj'\ 

and let v E he defined by vj = Tr{W*). 

The formula also depends on the choice of generators Ei, . . . ,6^^, where i' = n — 1, 
for the group U d ffp oi totally positive units. Using the regular representation 
p: U ^ r, defined via p{e) = W5{e)W~^, where 5{e) is the matrix diag(e'^^\ . . . , e^"-*), 
we identify the units Ej with elements Aj = pIejY G F. (Note that p is the row regular 
representation.) 

7.3. Using the bar notation 

[A,\---\A,]:= (1, A,, A,A2, . . . , ■ ■ ■ A,) G T", 

we have the following proposition expressing the zeta values in terms of the Eisenstein 
CO cycle: 

7.4. Proposition. Let Uf be the subgroup U (1 (1 + f), and let tt run through all 
permutations o/ {1, . . . , z/}. Then for s = 1, 2, 3, . . . , 

ab,f,l-s)=r] J2 $^sgn(7r)vl>([A.(i)|...|A.(,)])(p-\g,p(£)*t;). 

e&U/Ui 7T 

Here the sign rj = ±1 is determined by 

T] = {-ly sgn(det W) sgn(i?), 

where R = det(log£:p''), I < i,j < jy. 

Proof. This follows from Corollary, p. 595] by writing the fundamental cycle of 
in terms of the Aj . □ 

7.5. Example. We work out the above formula in the case of a real quadratic field 
F. Let £ > 1 be the fundamental unit of U, the group of totally positive units in F, 
and let Zwi + Zw72 = fb~^ be a Z-basis of fb~^. Such a basis determines a matrix 
A = (2 ^) G SL2{Z) and a vector v e Q'^ via 

/ EWi \fac\fwi\ 
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In addition, we get the normforms 

P{X) = N{xiWi + X2W2), Q{X) = N{xiW2 - X2W1). 
Let p be the smallest positive integer such that {A^ — l)v G Z^. Then 

(50) CF(b,f,l-s) = r7 Yl '^ihA){P''\Q,A'v), 

k mod p 

where rj = sgn(w2w|^^ — Wiw^2^), and, if s = 1, 

(51) d))(i'<2'"' = |;*("2' + |^2(ct.i-ot-2) 

(52) - J: S,il±^)mMl±^-v,), 

^-^ c c 

j mod c 

with an additional correction term — sgn(c)/4 on the right if 1; G 1? . The finite sum 
([5^ ) is the classical Dedekind-Rademacher sum (27ri)~^S' (Z^, ( J " ) , 1, t;) |^. Note that 

the number of terms in that sum equals |c| = |(e — — w')!, where w = W2/W1, 

and the prime is Galois conjugation. Depending on e, this number can be very large. 
To get a more efficient formula for calculating \E', we apply the euclidean algorithm 
to the first column of A an obtain a product decomposition 

A = B,---Bt, t>l, B,= ~J 

Then 

t-i 

(53) ^(1, A) = Y^iBi ■ ■ ■ Be)^{l, B^+i). 



£=0 

Here t is roughly log |c|. Thus the number of terms is effectively reduced from |c| to 
log |c|, since 



, 6 \ e^2s(t^_^p / 1 ^ A -^^2)e^^l+r2(v2) 



1 1 ; {2s)\ vol/ (l+n)! (l + rs)! 

where the rational numbers -Pr ( " (5 ) are the coefficients of the polynomial 

a p \ XI' x;' 



In the exceptional case s = 1, f = (1), the correction term 

— -{sgn(i(7 + h) + sgn(iy' + h)} 
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must be added to the right side of (|5^). 

As a numerical example, we choose F = Q(v^), e = (3 + v^)/2, f = b = (1), 
= -e, W2 = 1. Then r] = +1, A = V). ^ii(oi) = 3, while Psol??) = 2, 
-Pii ( 1 1 ) = Po2 ( 1 1 ) = 2. Hence, according to ( |50D and ([5^), we get for the value 
of the Dedekind zeta function of F at s = — 1, 

Cf{-1) = C((l), (1), -1) = 3(2 - 7 + 2)(-^) + 3(1)^ = 1. 

7.6. Example. As a second example, we consider the cubic field Q(^^) of discrimi- 
nant 148 given hj 6'^ — 6"^ — 36 + 1 = 0. According to [Q, the group of totally positive 
units U is generated by ei = -36"^ + 26* + 10 and £2 = 56*2 + 66' - 2. 

Let f = (2) and b = (1). Then fb'^ = Zwi + Zw2 + Zw^, where Wi = 2, W2 = 29, 
and W3 = 26"^. With respect to this basis, we find 



10 


3 


1 


2 


1 





-3 


-1 






Ai= { 2 1 ,^2 = 6 13 28 





Then v = (1/2, 0, 0)* and rj = 1. Let V C Q"^ be a complete set of representatives for 
the orbit of f + Z'^ under the action of U (via Ai and A2) on Q^/Z^. Note that V is 
a finite set. Thus 

(55) ab, f, 0) = 5^(^(1, A,, A,A2) - v|/(l, A2, A^A2)) (1, Q, v). 

Each term on the right of ( |5BD breaks up into 10 Dedekind sums: one of rank 
3, three of rank 2, and six of rank 1. Note that the rank 3 and rank 1 sums are 
diagonal, whereas the rank 2 sums are not. After making all sums diagonal, we find 
that each \E' in (^) has 30 terms. Applying the summation formula for diagonal sums 
(Proposition p.ll| ), we see that to evaluate \i/ on any element of V, we must sum 76 



terms. 

Since = 6162 = 1 mod f, we can take V = {v,A2v}. This yields 152 terms 
altogether, all of which sum to 0, and thus C((l)) (2),0) = 0. This agrees with 0, 
and also with the observation that since —1 preserves the congruence class of 1 + f, 
and the norm of —1 is —1, the special value at s = must vanish. 

Now let f = (3). Then we may take Ai,A2 as above, and v = (1/3,0,0)* and 
r] = 1. Since e^^ = eie^ = 1 mod f, we must sum 13 ■ 76 = 988 terms, and we find 
C((l), (3),0) = 2/3, again in agreement with 0. 

Note that to compute C((l); (A^)? 0) for various G Z, we must only compute 
Ai, A2, and thus "${1, Ai, A1A2) — \E'(1, y42, ^41^42) once. After this it is routine to 
compute special values at s = 0, and the complexity in (|55|) comes from ^V, which 
can be large, even for small values of A^. A table of C((l); (^)) 0) for several rational 
integers A^ is given below. The values for A^ = 2, 3, 5, 7 are also in H. 
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AT 


AT t 


AT 

i V 


AT t 


i V 


N t 
iV (, 


i V 


iV (, 


i V 


iV (, 






11 


-18 


21 


-78 


31 


74 


41 


382 


2 





12 


5 


22 


68 


32 


15 


42 


—228 


Q 
O 


o 
Z 


1 Q 

io 


— zz 


Zo 


1 o 
iZ 


QQ 
OO 


oz 


/I Q 


— ODD 


4 


1 


14 


-20 


24 


23 


34 


50 


44 


1 


5 


-4 


15 


42 


25 


106 


35 


-54 


45 


254 


6 


4 


16 


7 


26 


-24 


36 


-43 


46 


6 


7 


2 


17 


100 


27 


-190 


37 


20 


47 


-570 


8 


3 


18 


-32 


28 


25 


38 


22 


48 


-13 


9 


-10 


19 


82 


29 


242 


39 


156 


49 


-222 


10 


-2 


20 


4 


30 


6 


40 


2 


50 


178 



C((l), (A^),0) for the cubic field of discriminant 148. 



8. Witten's zeta function 

8.1. To give another illustration, we recall the definition of Witten's zeta function 
and show how our algorithms can be used to compute special values of this 

function at even integers. For unexplained notions from representation theory, the 
reader may consult [Q. 

Let g be a simple complex lie algebra, and let R be the associated root system. 
Let (respectively -R") be a subset of positive roots (resp. negative roots), and let 
A C be the set of simple roots. 

The roots R generate a lattice A/j in an ^-dimensional real vector space E endowed 
with an inner product ( , ). Let Ayj/ be the weight lattice, which is the dual of A^j 
with respect to this inner product. We denote by ^7 C -E the set of fundamental 
weights, which is the basis of Ky/ dual to A. 

Let W be the Weyl group of R. This is a finite group that acts on E via a refiection 
representation, and preserves the inner product and the lattices A/j and k.^^. There 
is a decomposition of E into a finite union of rational polyhedral cones, and W acts 
by permuting these cones. Let be the closed top-dimensional cone generated by 
^. 

Let n denote the set of isomorphism classes of complex irreducible representa- 
tions of g. It is known that elements of LI are in bijection with the set Ky^ fl C^, 
the dominant weights. Given A from this latter set, we denote the corresponding 
representation by tta. Then the definition of the zeta function associated to g is 

(56) Cg(5):= E (dimTrA)-^ 

AeAvKnc+ 

8.2. Let m > 1 be an integer. The special value C,^{^rn) can be computed using a 
Dedekind sum as follows. 
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Let p be one-half the sum of the positive roots. An apphcation of the Weyl char- 
acter formula Corollary 24.6] shows that for any dominant weight A, we have 



(57) dim7rA= J] 



(p + A,a) 



a&R+ ' 



It is known that any dominant weight A can be written as a nonnegative integral 
linear combination of the fundamental weights. Using this in (^), a computation 
shows that (|5^) becomes 



(58) Cs(2m)=M2™ 



ni-i(ai,a;)2™ 



Here M is the integer nae-R+(/^' c*^)) ''^ is the number of positive roots, and the G 
(Z>°)^ are the coefficients of the positive roots in terms of A. The pairing ( , ) is 
the usual scalar product on M^. 



8.3. We obtain a Dedekind sum by extending the sum (pSD to the whole lattice. 
Let C W be the span of the first i basis vectors, and let cr = (t(0) be an r x r 
integral matrix such that (crj,x) = (aj,x) for x G Z^, and such that deter = 1. Let 
e = {2k,..., 2k) G W. 



8.4. Proposition, (cf. [|, p. 507]; Cg(2fc) = fj5(Z^ ^(fl), e, 0). 



Hence these special values can be computed in polynomial time using our tech- 
niques. We conclude with two examples: sis and s^. We recommend verification of 
these formulas to the interested reader for a pleasant combinatorial exercise. 
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8.5. Proposition. Let ({s) be the Riemann zeta function. Then 

(59) ^a(2m)=8 2m -1 k«C(6m-.). 

0<j<2m ^ ^ 
j=0 mod 2 

(60) ^a(2m) = 16 5^ C"~!^0(^ + ^ + ^ + ^)' ^^^^^ 



0<i<2m 

(61) 



v-^ / 2m, + z — 7 — 1\ / 6m + « — t — 1\^, 

A = J]o<,<2m( ._/ ^^^^ k(j)C(t)C(12m-j-t), 



0<t<4m+j-j 
j,t=0 mod 2 



(62) 

B = i: „^,^.„ (""+,'7-')(''" + l-^-":\0KWC(12m-J-.), 

^ o<«<2m V yv4m + z- j- l/ 

j,M=0 mod 2 

(63) 

^ \ ^ /2m + i — A; — 1\ /6m + i — — f — 1\ ^ , ^ , , , 

^ = E 0,.,. ( C(^)C(.)C(12m-A.-.), 



0<D<4m+i— fe 
A:,i;=0 mod 2 



(64) 

E/2m + 2 — — 1\ /6m + 2 — — w — 1\ , , 



0<«;<2m 
k,w=0 mod 2 



8.6. Remark. The formula (|59|) was independently discovered by Zagier, S. Garoufa- 
lidis, and L. Weinstein p. 506]. 



8.7. Example. Here are some special values of Cats and (^h- 



2m 


(6m + 1)! ■ 6 ■ CB[3(2m)/(22™ • (27r)*^'") 


2 


1/(2-3) 


4 


19/(2- 3- 5) 


6 


1031/(3-7) 


8 


(ll-43-751)/(2-7) 


10 


(5-13-27739097)/(3-ll) 


12 


(17 - 29835840687589)/(3 - 5 - 7 - 13) 


14 


(2 - 17 - 19 - 89 - 127 - 6353243297)/7 


16 


(19 - 23 - 31 - 221137132669842886663)/(2 - 5^ - 13 - 17) 
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2m 


(12m + 1)! • (6m + 1) • (4m + 1) • 24 • CshC^m) / {W^"^ ■ (27r)^2'") 


2 


23/2 


4 


(3-7-14081)/2 


6 


(757409 • 23283173)/(5 ■ 7) 


8 


(3 • 11 • 1021 • 5529809 • 754075957)/2 


10 


(13 ■ 116763209 ■ 1872391681 ■ 3187203549787)/(5 ■ 11) 


12 


(17 • 1798397149 • 5509496891 • 6127205846988571484743)/ (3 •7-13) 
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